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ABSTRACT 

LSI +61° 303 is one of the few X-ray binaries with Be star companion from which both radio 
and high-energy gamma-ray emission have been observed. We present XMM-Newton and 
INTEGRAL observations which reveal variability of the X-ray spectral index of the system. 
The X-ray spectrum is hard (photon index F ^ 1.5) during the orbital phases of both high and 
low X-ray flux. However, the spectrum softens at the moment of transition from high to low 
X-ray state. The spectrum of the system in the hard X-ray band does not reveal the presence 
of a cut-off (or, at least a spectral break) at 10-60 keV energies, expected if the compact object 
is an accreting neutron star. The observed spectrum and spectral variability can be explained 
if the compact object in the system is a rotation powered pulsar. 

Key words: pulsars : individual: LSI +61° 303 - X-rays: binaries - X-rays: individual: 
LSI +61° 303 



1 INTRODUCTION 

The Be star binary LSI +61° 303 is one of the few X-ray binaries 
(along with PSR B 1259-63 and LS 5039) from which radio and 
very high-energy gamma-ray emission is observed. Radio emis- 
sion from the system is highly variable and shows periodicity of 
T = 26 .4960 days, which ca n be associated with the binary orbital 
period iOregorv et alj|20o3) . Radio observations reveal the pres- 
ence of 100 AU-scale jet in the syst em which places LSI +61 ° 303 
among the Galactic micro-quasars jMassi et alj[l993l l2004 . The 
system is also a Galactic "micro-blazar" du e to its association with 
100 MeV gamma-ray so urce 2CG 135+01 <Tavam et alll998h vis- 
ible up to TeV energies l' Albert et a l .120061) . 

iMendelso n & Maseh ( 1989) found an optical modulation in 
the V-band with the same periodicity as the radio periodicity. This 
is confirmed by Paredes et al. (1994); the optical modulation in 
the V-band is about 0.15 mag. Optical data allow to constrain the 
orbital param eters o f the system revealing the eccentricity of the 
orbit, e ^ 0.7 iCasares et al. 2 0o3) . However, the measurements are 
not sufficient to determine the nature of the compact object (neutron 
star or black hole), because the inclination of the orbit is poorly 
constrained. 

The X-ray emission from LSI +61° 303 is also variable and 
has nonthermal spectrum. Since the discovery, the system was 
twice monitored si multaneously in the X-ray and radio band s over a 
single orbital cycle iTavlor et alll996l : lHarrison et all2o'ool) . These 
simultaneous observations show that X-ray emission peaks almost 
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half an orbit before the radio. In both X-ray and radio bands the 
orbital phases of maxima of the flux "drift" from orbit to orbit. 

Two major types of models of radio-to-X-ray activity of 
LSI +61° 303 were proposed in the literature. Models of the 
first type, first introduced by iTavlor & Gregorvl il984h assume 
that activity of the source is powered by accretion onto the com- 
pact object. In the secon d class of models, first proposed by 
iMaraschi & Treve'3 il98lh . the activity of the source is explained 
by interactions of a young rotation powered pulsar with the wind 
from the companion Be star. 

If the system is an accreting neutron star or black hole, one 
expects to find a cut-off powerlaw spectrum in the hard X-ray 
band. The cut-off^ energy is normally at 10-60 keV iWhite et alj 
ll983t lFilippova etaL/OOa) for neutron stars and at ~ 100 keV for 
black holes I McClintock & Remillard 200^ (the cut -off can move 
to even higher energies in the "high/soft" state see e.g. lBelloni et alj 
i2006h ). Of course, it is possible that the emission from the "cen- 
tral engine" of the micro-blazar LSI +61° 303 is "masked" by the 
emission from the jet, which can be beamed toward an observer 
on Earth, similarly to what is observed in blazars. In this case the 
hard X-ray spectrum of the source is a superposition of the pow- 
erlaw jet emission and emission from the accretion disk. If the jet 
and accretion contributions to the X-ray spectrum are comparable, 
then emission from the accretion disk should at least produce an 
observable spectral feature (e.g. a bump, a break or turnover) in the 
10-100 keV energy band. 

Below we present a study of LSI +61° 303 in the 0.5-100 keV 
energy band with XMM-Newton and INTEGRAL which shows 
that the spectrum of the source is well fit by a simple powerlaw, 
without any signature of high-energy cut-off. Moreover, the X- 
ray powerlaw matches smoothly to the higher-energy spectrum in 
the 100 keV - 10 MeV band found with the CGRO instruments 
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OSSE jTavanietalj|l996h and COMPTEL Jvan Diik et ai]|l99el) . 

The broad band spectram of the source does not exhibit a cut off 
up to at least GeV energies. The spectrum of LSI +61° 303 in the 
X-ray - soft gamma-ray band is thus very different from the typical 
spectra of accreting neutron stars and black holes. 

Contrary to accreting compact object models, a featureless 
powerlaw keV - MeV spectrum is expected in the "rotation pow- 
ered pulsar" scenario. Unfortunately, the population of X-ray bina- 
ries with a rotation powered pulsars as compact objects is not so 
rich as that of accreting pulsars. In fact, the only firmly established 
example of X-ray binary with a yo ung radi o pulsar as a compact 
object, PSR B 1259-63 Johnston et aljl992h . exhibits a featureless 
powerlaw spectrum in the X-ray band which smoothly continues 
to the higher energies, exactly as in LSI +61° 303. Moreover, the 
overall spectral energy distributions and variability properties of 
the two systems are qualitatively very similar: large radio outbursts 
which occur once per orbit, similar X-ray spectra and high-energy 
gamma-ray emission extending up to the TeV energy band. Recent 
observations of PSR B 12 59-63 during its 2004 periastron passage 
iChernvakova et al]|200(il) show that the X-ray emission from the 
system is most probably produced via inverse Compton scattering 
of the UV radiation from the Be sta r by the pulsar wind electr ons, 
the same mechanism as proposed bv lMaraschi & Treved Jl98lh for 
LSI +61° 303 

The XMM-Newton and INTEGRAL observations of 
LSI +61° 303 reported here, reveal the variability of the X- 
ray spectrum of the system during the transition from high to low 
X-ray state. We find that the spectra of the system in high and low 
X-ray states are very similar, but during the transition to the low 
flux state the spectrum softens. It turns out that explanation of such 
behavior of the system is not trivial and, in fact, the shape and 
variability of the spectrum put tight constraints on the injection 
mechanism of high-energy electrons responsible for the X-ray 
emission as well as on the geometry of the X-ray emission region. 

This paper is organized as follows: in Section 2 we describe 
the details of the INTEGRAL and XMM-Newton data analysis. The 
results (imaging, spectral and timing analysis) are presented in Sec- 
tion 3. In Section 4 we discuss the physical implications of the re- 
sults for the synchrotron - inverse Compton broad band model of 
the source. This gives us the possibility to constrain physical pa- 
rameters of the system and develop a physical model of the source 
which we study in more details in Section 5. Finally, we summarize 
our findings in Section 6. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 INTEGRAL observations 

Since the launch of INTEGRAL <WinkIer et al]|2003h on October 
17, 2002, LSI +61° 303 was several times in the field of view of 
the main instruments during the routine Galactic plane scans and 
pointed observations (see Table Q for details). Most of the times 
the distance of the source from the center of the field of view was 
too large to use the X-ray monitor JEM-X. Therefore IBIS/ISGRI 
iLebrun et all2003h is the only instrument we can use in our anal- 
ysis of this source. In this analysis we have used the version 5.1 
of the Offline S cience Analysis (OSA) software distributed by the 
ISDC <Courvoisi er et aL2Q03). 

In our analysis we have used all available public data spread 
over the period from the January 2003 (rev 25) until March 2005 
(rev 288). Overall we have analyzed 600 science windows which 
resulted in an effective vignetting corrected exposure of 273 ksec. 



Table 1. Journal of the INTEGRAL obsei-vations of LSI +61° 303. 



Data 


Orbital phase 


Effective 


20-60 keV Flux 


r 


Set 




Exposure (ks) 


10-"ergs s-'cm-^ 




11 


0.4 - 0.6 


50 


3.8+0.6 


1.7+0.4 


12 


0.6-0.8 


23 


3.0+1.0 


3 6+'-'' 


13 


0.8-0.4 


200 


2.4+0.3 


1.4+0.3 


lav 




273 


2.5+0.3 


1.6+0.2 



■ Given errors represent 68% confidence interval uncertainties. 



Table 2. Journal of XMM-Newton observations of LSI +61° 303 



Data Observational 


Date 


MJD 


Orbital 


Exposure (ks) 


Set 


ID 




(days) 


Phase MOS1/MOS2/PN 


XI 


0112430101 


2002-02-05 


52310 


0.55 


5.8/5.8/5.0 


X2 


0112430102 


2002-02-10 


52315 


0.76 


5.8/5.8/5.0 


X3 


0112430103 


2002-02-17 


52322 


0.01 


1.5/1.5/5.0 


X4 


0112430201 


2002-02-21 


52326 


0.18 


7.0/7.0/5.1 


X5 


0112430401 


2002-09-16 


52533 


0.97 


6.2/6.2/6.0 



In order to study the spectral variations of the source we have 
grouped the data into three parts covering orbital phases 0.4 - 0.6, 
0.6 - 0.8 and 0.8 - 0.4 (see Table0. The spectral analysis was done 
with the use of the mosaic_spec tool, which extracts spectra for a 
given sky position from mosaic sky images. 



2.2 XMM-Newton observations 

XMM-Newton has observed LSI +61° 303 with the EPIC instru- 
ments five times during 2002. Four observations have been done 
during the same orbital cycle, and the fifth one has been done seven 
months later. The log of the XMM-Newton data analyzed in this 
paper is presented in Table|2| 

The XMM-Newton Observation Data Files (ODFs) were ob- 
tained from the on-line Science Archive' ; the data were then pro- 
cessed and the event-lists filtered using xmmselect within the Sci- 
ence Analysis Software (sas) v6.0. 1. In all observations the source 
was observed with the PN detector in the Small Window Mode, 
and with MOSI and M0S2 detectors in the Full Frame Mode. In 
all observations a medium filter was used. 

The event lists for spectral analysis were extracted from a 40" 
radius circle at the source position for MOSI, from a 30.0" radius 
circle for MOS2, and for PN a region of 25" around the source 
position was chosen. Background photons were collected from a 
region located in the vicinity of the source with the same size area 
as the one chosen for the source photons. 

For the spectral analysis, periods affected by soft proton flares 
need to be filtered out. To exclude them we have used script 
xininlight.clean. csh The source countrate in 2 - 10 keV en- 
ergy range varied from about 0.03 to 0.06 cts/s for MOS instru- 
ments, and from 0.05 to 0.10 cts/s for PN. With the help of the 
GRPPHA FTOOL we have re-binned the spectra, so that each energy 
bin has at least 30 counts. 

Data from MOSI, M0S2 and PN detectors were combined in 
the spectral analysis to achieve better statistics. 



http://xmm.vilspa.esa.es/extenial/xmm_data_acc/xsa/index.shtml 
http://www.sr.bham.ac.uk/xmm2/xmmlight_clean.csh 
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Figure 1. Contour plot of the XMM-Newton field of view for the XI obser- 
vation. A total of 5 contours were used with a square root scale between 1 
counts per pixel (outer contour) and 25 counts per pixel. 
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Figure 2. Contour plot of the INTEGRAL significance map of the 
LSI +61° 303 region in 22.1-63.2 keV energy range with a 273 ks effec- 
tive exposure. A total of 6 contours were used with a linear scale between 3 
o" (outer contour) and 8 cr (inner contour). 



3 RESULTS 

3.1 Imaging Analysis 

Three sources were clearly detected in the XMM-Newton field of 
view during its LSI -1-61° 303 observation program (see Fig.Qfor 
the contour plot of XMM-Newton field of view for the MJD 52310 
(XI) observation). Besides LSI -1-61° 303 itself we have detected 
two more sources that were also detected by ROSAT, and were 
identified as two optical sources BD -1-60 536, and HD -1-16 429 
bylOpldoni & Mereghetti tl99l) . In the X-rays BD -h6 536 has 
also been detected by ASCA in 1994 <Leahv et alJll997h . HD -1-16 
429 has also been observed in radio during d eep VLA observations 
of the LSI -f61° 303 field (Marti et a lll998l) . 

The presence of a quasar QSO 0241-1-622 some 1.4° away 



from LSI -1-61° 303 was a problem for the previous (non-imaging) 
observations of LSI -1-61° 3 03 above 20 keV. Care ful analysis of 
RXTE/REXTE data allowed [Harrison et all hoOOh to distinguish 
the two sources. The imaging capabilities of the INTEGRAL (sev- 
eral arcminutes angular resolution) allow, for the first time, to 
clearly separate the two sources. Both sources are clearly seen in 
the 22.1-63.2 keV INTEGRAL mosaic of the field (Fig|2j, at 8.1 cr 
(LSI -t-61° 303), and 12.7 a (QSO 0241-1-622) detection level. QSO 
0241-1-622 is several times brighter than LSI -1-61° 303. BD -1-60 
536, and HD -h16 429 detected by XMM-Newton are too faint to be 
observed with INTEGRAL. 

3.2 Timing analysis 

Most previous papers on LSI -1-61° 303 ha ve defined the ra- 
dio phase according to the ephemeris of either ll^vlor & Gregory! 
( 1984), or Gregorv et al. ( 1999). In the cur rent paper we use the 
latest ephemeris of the system derived by iGregorv et alj i2002h 
(Port = 26.4960 d, ro=JD 2,443,366.775). This allows a consis- 
tent treatment of current and historica l observations. Optical ob- 
servations allowed lCasares et alj i2005h to build an orbital solution 
of the LSI -1-61° 303 system. This solution constrains the perias- 
tron phase to (po = 0.23 ± 0.02, and the eccentricity of the system 
e = 0.72 ±0.15. 

Third, forth and fifth (from the top) panels of Fig. |^ show 
the dependence of the 20-60 keV, 2-10 keV and 0.5-2 keV X-ray 
fluxes on the orbital phase. Besides the new INTEGRAL and XMM- 
Newton dat a we show in this Figure all availabl e archival data 
from ASCA jLeahv et ajl997h . RXTE and ROSAT jHarrison et all 
2000). 

To allow straightforward comparison with the data in other 
energy bands we also show in Fig.|3]the lightcurves of the source 
in the broad (radio-to-gamma-ray) energy interval. The upper- 
most panel shows the highest energy (TeV) lightcurve reported in 
Albert et al, 1 2006). The second panel shows the EGRET (GeV) 
data (ITavani et alj|l998h . Finally, the bottom panel of Figure |^ 
shows radio data sim ultaneous to RXTE and ROSAT observations 
jHarrison et aj200(j) . 

In agreement with previous observations the 2-10 keV flux 
during XMM-Newton observations was observed to be variable with 
phase. It was rather low at zero phase, reached its minimum at pe- 
riastron (near phase 0.2), and became approximately three times 
higher at phases larger than 0.5. Comparison of the X-ray and ra- 
dio lightcurves of different years shows that the X-ray flux from 
the system varies in a more regular way than the radio flux: the X- 
ray data points from different orbital cycles lie more-or-less on the 
same curve, with the only exception of the second A5CA point. 

3.3 Spectral Analysis 

The spectral analysis was done with NASA/GSFC XSPEC soft- 
ware package. In Fig.|31the XMM-Newton spectra for the XI, X3 
and X4 observations , as well as the average INTEGRAL spectrum 
of the source are shown. We do not show the spectra of observations 
X2 and X5 because they are very similar to that of XI. 

For comparison we also give in this Figure the 20-100 keV 
spectrum of QSO 0241-t-622 (the nearby source which was difficult 
to distinguish from LSI -t-61° 303 in earlier non-imaging observa- 
tions). The spectrum of QSO 0241-1-622 is well fitted by a power- 
law with the photon index Tgso = 2. 1 ± 0.2 and 20-100 keV flux 
F = (5.5 ± 0.4) X 10"" erg/(cm^s). The source did not show any 
variability during the 2-year period of INTEGRAL observations. 
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Figure 3. Comparison between the TeV (top), GeV gamma-ray orbital lightcurves of LSI +61° 303 with the hard X-ray (20-60 keV), X-ray (2-10 keV), soft 
X-ray (0.5-2 keV) and radio (bottom) orbital lightcurves. The 0.5 - 2 keV and 2-10 keV X-ray flux is given in 10"" ergs cm^^s"' . GeV y-ray flux (> lOOMeV) 
is given in the units of lO^^ph cm^^s"' . TeV flux (> 400 GeV) is in the units of 10"" ph cm"^s"'). To guide the eye we have connected with lines data from 
the same orbital cycle. The time is assumed to increase from left to right, thus if the observation has started at the end of the orbital cycle it is shown with two 
lines of the same color, the line starting at larger orbital phase coimects data points taken earlier 
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Figure 4. XMM-Newton PN LSI +61° 303 spectra from the XI, X3 and 
X4 observations. For INTEGRAL the averaged 2003 - 2004 spectrum is 
shown. We also show on this plot the spectium of the neai'by source QSO 
0241+622. 



A simple power law with photoelectric absorption describes 
the spectrum of LSI +61° 303 well, with no evidence for any line 
features. In Table |3] we present the results of the three parameter 
fits to the XMM-Newton data in the 0.5 - 10 keV energy range. The 
uncertainties are given at the Icr statistical level and do not include 
systematic uncertainties. The values for the flux and the powerlaw 
index were obtained assuming that the hydrogen column density 
Nn was the same for all the datasets. The value = (0.49 ± 
0.02) X 10"" cm"^ was found from the simultaneous fit of all the 5 
XMM-Newton datasets. 

The graphical representation of the evolution of the spectral 
parameters along the orbit is given in Figure|5|which includes soft 
X-ray from XMM-Newton and ASCA observation s (top panel), and 
averaged radio data obtained bv lRavetalJjl997h with the National 
Radio Astronomy Observatory Green Bank Interferometer (GBI) 
(bottom panel). 

In the XMM-Newton energy band the powerlaw photon index 
is r 1.5 in all the observations except for the observation X3 
where we find significantly softer spectrum with T = 1.78 + 0.04. 



INTEGRAL and XMM-Newton observations of LSI +61° 303 5 



1.6 



1.4 



1.2 - 



YMM 0Z.200Z 
Xm 09.2002 
ASCA 02.1994- 




0.5 1 

Orbital phase 

Figure 5. Evolution of the photon index F as a function of the orbital phase 
for XMM-Newton and ASCA (top panel), and radio GBl (bottom panel) ob- 
servations of the LSI +61° 303 system. 



Table 3. Spectral parameters for XMM-Newton Observations of 
LSI +61° 303.* 



Data 


F(2-10keV) 


r 


Nh 


(dof) 


Set 


10""erg s"' 




{ cm-2) 




XI 


1.30±0.04 


1.60+0.03 


0.49+0.02 


474(456) 


X2 


1.24±0.04 


1.54+0.03 


0.49+0.02 


440(446) 


X3 


0.61+0.03 


1.78+0.04 


0.49+0.02 


184(188) 


X4 


0.44+0.03 


1.52+0.06 


0.49+0.02 


137(150) 


X5 


1.25+0.05 


1.55+0.02 


0.49+0.02 


488(496) 



* Given eiTors represent 68% confidence interval uncertainties. 

This indicates that the X-ray spectrum softens during the transition 
from higher to lower flux state. 



4 BROAD BAND SPECTRUM OF THE SOURCE. 

The rich multiwavelength variability properties of the system, pre- 
sented in Figs.|3||4||5lsuggest that the physical model of the system 
which is able to explain all the observed spectral variations should 
be very restrictive. However, up to now most of the model building 
for the source concentrated on explaining the details of either ra- 
dio or X-ray or gamma-ray variability, paying little attention to the 
broad-band properties of the source. 

As mentioned in the introduction, historically there are two 
broad classes of physical models of the source activity. The first 
one assumes that the activity is powered by accretion onto the com- 
pact object (either neutron star or a black hole) (Tay lor & Gr egor^ 
Il984h . Accreting X-ray binaries with a Be star companion con- 
stitute the majori ty of the Galactic high-mass X-ray binaries 
iZiolkowskill2002h . Most of the Be star X-ray binaries contain an 
accreting neutron star as the compact object and most of them are 



transient sources. The spectra of these sources are characterized 
by the presence of an e xponential cut-off i n the hard X-ray ban d 
at energies 10-60 keV iWhite et alj|l983l: iFilippova et ai]l2005h . 
Our analysis of the INTEGRAL data shows that in the case of 
LSI -1-61° 303 no high-energy cut-off is found at the energies be- 
low 100 keV, the fact that does not fit well into the "conventional" 
accretion scenario. It is possible that in the case of LSI +61° 303 
emission from the accretion disk is "masked" by the emission from 
the jet, which can have a simple powerlaw spectrum. If accretion 
flow and jet would give comparable contributions in the hard X- 
ray band, one would still expect to observe a spectral feature be- 
tween 10 and 100 keV which would result from the disappearance 
of the accretion component at higher energies. The non-observation 
of such feature in our data indicates that the luminosity of the accre- 
tion flow in LSI -1-61° 303 is much lower than the observed X-ray 
luminosity, Laccction ^ lO^^"* erg/s. Since, by assumption, the source 
activity (including the jet) is powered by accretion, the large ap- 
parent luminosity of the jet (Lj^, ~ 10'^ erg/s, of the order of the 
gamma-ray luminosity) can be explained by a large beaming factor 
of the jet, 6-2. 

Radio observations suggest that the jet is moderately relativis- 
tic, with the speed /3 ~ 0.4 0.6. It is presumably precessing 
and the angle with respect to the line of sight v aries from 9^0° 
jMassi et all200lh to 6) ^ 78° jMassi et all2004h . This implies that 
the Doppler factor of the jet, 6 = -^/l -fi^/ {l±f3 cos 9), varies in the 
range I < 6 < 2 along the orbit. If the high gamma-ray luminosity 
of the source during a fraction of the orbit is explained by the large 
Doppler factor of the jet (Lj>, ~ 6'*), one would expect that the jet 
contribution drops during the orbital phase with the low Doppler 
factor. In this case the "accretion" contribution would show up in 
the low X-ray flux state. However, our data show that the spectrum 
of the source is well fit by a featureless powerlaw with the same 
photon index both in the high and low X-ray states. It is not impos- 
sible that the accretion component is, similarly to jet component, 
modelled by a powerlaw, but in this case the coincidence of the 
photon indexes in the high and low flux states implies an artificial 
"fine tuning" of parameters. Although we find little support for the 
accretion model in our data, further multiwavelength observations 
are necessary to better constrain the Doppler factor and the hard 
X-ray spectral properties of the source low state. 

The timing and spectral properties in the keV-MeV energy 
band favor the second class of models, which assume that the sys- 
tem activity i s powered by the rotation energy of a relatively young 
radio pulsar iMaraschi & Treve3ll98lh . At least one more Be star 
X-ray binary which is powered in this way is known, namely, 
PSR B1259-63. In fact, the spectral properties of PSR B1259-53 
system are very similar to the LSI +61° 303: strong radio flares 
and non-thermal X-ray and gamma-ray emission are detected dur- 
ing the periods of pu l sar passing through the di s k of the Be star 
ilohnston et alj 120051: IChernvakova et alj 120061: I Aharonian et alj 
l2005t) . In the model of iMaraschi & Treve J jl98ll) the radio emis- 
sion is the synchrotron emission from the pulsar wind electrons 
while the X-rays are produced via the inverse Compton scatter- 
ing of the soft photons from the Be star by the same electrons. 
Roughly the same conclusion was dedu ced from the recent XMM - 
Newton observations for PSR B 1259-63 jChemvakova et alj2006l) . 
The main differences between the two systems are that (a) the pul- 
sar itself was never detected in the case of LSI +61° 303 and (b) in 
the case of PSR B 1259-63 no AU-scale jets were seen (one should 
note that PSR B 1259-63 is in the Southern hemisphere where no 
VLBI observations are available). 

In the case of the "young pulsar" scenario, radio. X-ray and 
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Figure 6. Broad band spectrum of LSI +61° 303 . Radio data points are taken from lStrickman et alj 1 1993) . CGRO data points are from lTavani et alj il99d 
[l998). MAGIC data points are from lAlbert et al] l20Qd) . The solid (dashed) line shows the model fit within the synchrotron-inverse Compton model for the 
high (low) flux state of the source. The dotted line shows possible contribution from the proton proton interactions. The values of parameters of the model fits 
are cited in the text. 



gamma-ray emission come from the region where the pulsar wind 
interacts with the wind of Be star. The high-energy particles re- 
sponsible for the synchrotron and inverse Compton emission can 
either originate from the relativistic pulsar wind or be accelerated 
in the shock formed at the interface of the stellar and pulsar wind. 

In the inverse Compton model of X-ray emission the soft pho- 
tons upscattered to the X-ray energies are the UV photons emitted 
from the Be star. If the distance from the X-ray emission region to 
the star is R, the energy density of UV photons at the location of 
the pulsar can be estimated as 



4nR^c 



2 X 10^ 



IC* erg/s 



10'^ cm 



R 



erg/cm'' (1) 



where L, is the luminosity of the Be star. Relativistic electrons up- 
scatter the UV photons up to energies 



keV 



(2) 



L2x 10^ KJLlOMeV 

where T is the temperature of the Be star Apart from producing 
inverse Compton emission in the X-ray band, the 1-10 MeV elec- 
trons emit radio synchrotron radiation. The typical frequency of the 
synchrotron photons is 



7 



B 

Tg 



GHz 



(3) 



10 MeV J 

where B is the magnetic field strength in the emission region. In 
the synchrotron - inverse Compton model the ratio of the radio and 
X-ray luminosities of the system, L,-adio/ix-ray, is equal to the ratio 
of the energy densities of magnetic field, 

(/s-4x 10-'[^] erg/cm' 
to the energy density of radiation Q: 

^radio 



(4) 



Ub 

Urad 



2 X 10" 



B 

Tg 



R 



W- cm 



lO^** erg/s 



(5) 



The last relation imposes a restriction on the possible values of 



magnetic field in the emission region. Taking into account that the 
distance from Be star to the compact companion is constrained by 
the binary orbit parameters roughly in the limits lO''^ cm < < 
10" cm and taking the observed ratio Liadio/^'X-ray ~ 10"', the 
magnetic field is constrained to be 

B < {several) G (6) 

An example of synchrotron - inverse Compton model fit of the 
broad band (radio to gamma-ray) spectrum of the source is shown 
in Fig.|6| The model parameters used for the fits of the broad band 
spectrum of the system in high (low) flux state in Fig. |6| are 6 = 
0.35 G (S = 0.25 G), the electron spectrum described by a broken 
powerlaw with the spectral index = 2 below the break energy 
Etr = 2 X lO** MeV (£,„. = 10** eV) and T, = 3.5 (T, = 4) above 
the break energy. In the case of low radio state the break is at a bit 
lower energy and the electron spectrum is softer above the break 
energy, to fit the radio spectrum of the source which clearly exhibits 
a cut-off at ~ 10 GHz frequencies. One can see that in this case the 
inverse Compton spectrum softens at the energies £ ~ 100 keV. 

It is interesting to note that even if the electron spectrum 
extends to energies higher than 10 GeV, one still expects to ob- 
serve a softening of the gamma-ray spectrum at 10 GeV energy 
because above this energy the inverse Compton scattering proceeds 
in the Klein-Nishina regime and if the electron spectrum is a pow- 
erlaw with the spectral index Tg, the inverse Compton spectrum 
below 10 GeV is a powerlaw with the photon index (F^ -I- l)/2 
(dNy/clE ~ £:-('"«+i)/2) while above 10 GeV the spectrum softend 
to roughly clNy/dE ~ j?""^'^'* I ni?. The spectrum above 400 GeV 
measured by MAGIC telescope jAlbert et al J200d) is indeed softer 
than the EGRET spectrum in the GeV band. This can be an indi- 
cation of transition to the Klein-Nishina regime. However, the ob- 
served softening is less than expected in this case and it is possible 
that emission above 100 GeV is dominated by anoth er mechanism 
(e.g. b y pion decay gamma-rays, as suggested by iRomero et alj 
i2005h ). We come back to this possibility in Section 5. 

The X-ray luminosity of the system is several orders of mag- 
nitude higher than radio luminosity. This enables to conclude that 
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the inverse Compton energy loss always dominates over the syn- 
chrotron energy loss. Indeed, the inverse Compton cooling time is 



t,c - 10" 



10" erg/cm 



U,r, 



1 keV 



1/2 



while the synchrotron cooling time is 



^svnch ~ 10 



1 G 


3/2 


1 GHz 


1/2 
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B 









(7) 



(8) 



Factor of 10 variations of the inverse Compton (X-ray) lumi- 
nosity of the source can be produced if the typical distance from the 
emission region to Be star varies by a factor of 3. This is smaller 
than the variations of the distance from the compact object to Be 
star, estimated from the eccentricity of the binary orbit. But, in fact, 
in the models of interaction of the pulsar wind with the stellar wind 
from the companion star the bulk of emission from the system is 
supposed to come from the shock region at the contact surface of 
the two winds, which is situated in between the star and the com- 
pact object and possibly extends to the distances larger than the 
binary separation. 

In order to explain the observed spectral properties of the 
source, one has to compare the inverse Compton cooling time to 
the escape time from the emission region. Several velocity scales 
are present at the interface of stellar and pulsar winds. If the pul- 
sar wind mixes with the stellar wind, the escape of the high-energy 
particles could be slowed down to the velocity of the "slow" equa- 
torial wind from Be star, vvu^. ~ 10* -r 10^ cm/s. The same esti- 
mate of the escape velocity is obtained if the high-energy parti- 
cles are stellar wind particles accelerated in the shock at the in- 
terface with the pulsar wind. To the contrary, if the two winds 
flow along the contact surface without significant mixing, the es- 
cape velocity of the pulsar wind particles is comparable to the 
speed of light, Vp„/„r ~ 10'° cm/s. Correspondingly, two differ- 
ent scenarios are possible: that of a "slow escape", with the es- 
cape time fj,c ~ R/vjjsi, ~ 10* s and that of a "fast escape", with 

^esc ^ I pulsar ~ 10 S. 

In the case of "slow escape", the inverse Compton emission 
should be produced by the cooled electron population. For a pow- 
erlaw distribution of electrons, the steady state spectrum has a pho- 
ton index F^ which is softer than the spectral index Fjj„j of injection 
spectrum, F^, = F(, i„j -I- I. The observed photon index in X-ray band 
is Tpi, = (Fp + I)/2 1.5. This implies that Fp i„j -1-1^2, which 
means that the injection spectrum should be very hard, F^ inj - 1. 
Such spectrum hardly can be produced by shock acceleration at 
the contact surface of pulsar and stellar winds. Thus, in the "slow 
escape" scenario the only possibility is that electrons are initially 
injected only at high energies (that is, they originate from the cold 
pulsar wind with large bulk Lorentz factor) and the inverse Comp- 
ton cooling would leads to the formation of "universal" cooling 
spectrum with F^ = 2, regardless of the details of the injection 
spectrum at high energies. 

In the case of "fast escape" the high-energy electrons also 
originate in the pulsar wind (by assumption). In this case the in- 
verse Compton cooling time at 1-100 keV energies is comparable 
to the escape time and, in the case of powerlaw injection spectrum, 
one expects to observe a "cooling break" (softening by AF,,/, 0.5) 
at the energy at which the escape time is equal to the cooling time. 
From the model fit to the broad-band spectrum, shown in Fig|6| one 
can see that a break could be situated at the energies above MeV 
(electron energies £f j,- ~ 100 MeV) which is, in principle, com- 
patible with the "fast escape" scenario. In fact, the "fast escape" 



scenario is not so different from the "slow escape" one, because the 
spectrum of the high-energy electrons from the pulsar wind with 
large Lorentz factor has a low-energy cut-off. Escaping electrons 
cool down only to the characteristic energy E^jir ~ 100 MeV in the 
direct vicinity of the pulsar and Be star. As a result, 100 MeV elec- 
trons are injected in a larger region around the binary system and 
continue to cool (at a longer time scale) forming a characteristic 
E^- powerlaw spectrum below 100 MeV. 

The observations presented in the previous sections clearly re- 
veal the softening of the spectrum during the decrease of the X-ray 
luminosity. The softening of the spectrum can be explained assum- 
ing that the injection of electrons at high energies drops and higher 
energy electrons (emitting in the hard X-ray band) cool faster than 
the low energy ones (emitting in the soft X-ray band). Note that 
this explanation implies that the size of the region from which the 
bulk of X-ray emission comes is at least R ~ 10'^ cm, otherwise 
the inverse Compton cooling time would be too short to explain 
the softening on day time scale. 

Although radio and X-ray emission are produced by the same 
population of electrons, they are not necessarily produced in the 
same spatial region. The strongest synchrotron emission is pro- 
duced in the region with higher magnetic field which is situated in 
the direct vicinity of the pulsar, while the strongest inverse Comp- 
ton emission is produced in the region with highest soft photon 
background density, which is situated close to the companion star. 
This means that, in general, one does not expect a strong corre- 
lation of the X-ray and radio lightcurves. Taking into account the 
fact that 10-100 MeV electrons can propagate over a significant 
distance during the inverse Compton cooling time 0, the detailed 
modelling of the X-ray and radio lightcurves requires numerical 
simulations, which would take into account (somewhat uncertain 
and time-variable) geometrical properties of the system, such as, 
e.g., the geometry and density profile of the equatorial disk of Be 
star and its inclination to the orbital plane. We leave this for the 
future work. 



5 PHYSICAL MODEL OF THE SOURCE. 

The above modeling of the broad-band spectrum of the system 
shows that the physical model of the source (which we assume to 
be based on the "young pulsar powered source" model, see above) 
should explain several key features derived from the analysis of ra- 
dio. X-ray and gamma-ray data. 

In particular, we have seen that the X-ray spectrum and spec- 
tral variability are satisfactory explained if the relativistic electrons 
responsible for the X-ray emission are initially injected at high en- 
ergies and subsequently cool down forming the characteristic E^^ 
cooling spectrum. The physical model should provide an explana- 
tion for the electron injection only at high energies. The analysis of 
the radio and X-ray variability of the system shows that the injec- 
tion of high-energy electrons should be variable. 

Injection from high energies can be explained if the electrons 
responsible for the X-ray emission originate from the cold pulsar 
wind with bulk Lorentz factor > {several) x 10^. In this case all 
electrons have initial energies larger than 100 MeV and the electron 
spectrum below 100 MeV formed in the process of inverse Comp- 
ton cooling has the spectral index F^ ^ 2. However, an immediate 
difficulty with such simple injection model is that the injection rate 
is not expected to vary with time, contrary to what is observed. 

The variable injection rate could be provided if relativistic pro- 
tons are either present in the pulsar wind or accelerated in the shock 
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at the contact surface of pulsar and stellar wind. Such relativistic 
protons could interact with the low energy protons from the disk 
and produce injection of electrons at energies above ~ 100 MeV 
via production and subsequent decays of charged pions. The pion 



production cross section is CTp 



10-2* cm^. The life time of 



high-energy protons in the disk with the density n ~ 10 cm" 
is tpp = {/(cTppnc) ^ 10* s, comparable with the escape time from 
the disk. This means that (a) protons can efficiently transmit their 
energy to the products of pion decay (gamma-rays, neutrinos, elec- 
trons) during their propagation through the disk, (b) proton-proton 
interactions can provide the source for injection of 10-100 MeV 
electrons over the large AU-size region, from which the bulk of X- 
ray emission supposedly comes and (c) variable injection rate of 
high-energy electrons is explained by the variations in the density 
of the stellar wind protons along the pulsar orbit. 

A straightforward consequence of high-energy proton interac- 
tions in the Be star disk is the appearance of additional component 
in the high-energy y-ray spectrum, resulting from the two photon 
decays of neutral pions. The spectrum of the pion decay gamma- 
rays in the GeV-TeV energy band has the same spectral index as 
the spectrum of the high-energy protons. In Fig.|6|we show possi- 
ble contribution to the y-ray spectrum of the source which can be 
produced in result of proton-proton interactions. 

If the system activity is powered by a young pulsar, the pulsar 
should be also responsible for the formation of the 100 AU scale 
jet in the system. The fact that rotation powered pulsars can pro- 
duce jets is beyond doubt, since such jets are observed in the pul- 
sar wind nebulae of e.g. Crab and Vela pulsars. In the case of the 
pulsar interacting with the wind of companion star a formation of 
even more complicated ex pend ed emission features which can im- 
mitate the jets is possible ( Dubus 200(3'). However, no "immitation" 
of the jet is really needed, because the mechanism of formation of 
the jet in the system of a young pulsar interacting with a stellar 
wind can be the same as in the system of a young pulsar interacting 
with interstellar medium. The naive qualitative estimate shows that 
anisotropic power injection into the stellar wind should, in general, 
lead to the formation of an outflow from the system with param- 
eters close to the expected physical parameters of the jet. Indeed, 
the total energy which is injected during the escape time scale is 



10 erg. Such kinetic energy transmitted to a vol- 



ume of the size of about the thickness of the disk, H 
with the density n ~ 10'" cm"'' (total mass is M 
pushes the matter in the volume to move with the speed 



10'" cr 
nipjiH^ -10-^ 



O.lc 



1035 erg/s 



10* s 



10" 



10'- cm 



H 



(9) 

This velocity is of the order of the velocity of the fast polar wind 
from Be star and is much larger than the typical velocity inside the 
disk, Vrf,5t ~ 10* cm/s. Qualitatively, it is possible that (anisotropic) 
injection of large kinetic energy of the pulsar wind into the disk 
is responsible for the formation of a jet-like outflow (or "expand- 
ing bubbles" implied in the models of radio flares) from the disk 
in which relativistic protons and electrons from the pulsar wind 
are mixed with the low-energy protons from the disk. The proton- 
proton interactions resulting in production of pions which subse- 
quently decay and inject high-energy electrons can continue also 
in the outflow thus providing a mechanism of injection of high- 
energy electrons at larger scales up to the scale of the ~ 100 AU jet 
observed in the system. 



6 SUMMARY. 

In this paper we have presented the INTEGRAL and XMM-Newton 
observations of LSI +61° 303. We have found that the overall spec- 
trum of the system in 0.5-100 keV band is well fit by a featureless 
powerlaw with the photon index Tpi, ^1.5 both in the high and low 
flux states. The 0.5-100 keV powerlaw spectrum matches smoothly 
the higher-energy spectrum of the source in 100 keV - GeV band. 
The powerlaw with the same photon index seems to continue with- 
out a cut-off up to 1-10 MeV energies. Non-observation of a cut-off 
or a break in the spectrum at 10-100 keV energies, typical for the 
accreting neutron stars and black holes, implies that the accretion in 
the source proceeds at a very low (if not zero) rate. This favors the 
scenario in which the compact object is a rotation-powered pulsar. 
We have discovered that the spectrum is hard in both high and low 
X-ray states, but softens during the transition from the high flux to 
the low flux state. The spectral characteristics of the system favor 
the model in which electrons are initially injected at high energies 
and subsequently cool to lower energies forming the characteris- 
tic "cooling" spectrum with the spectral index 2. We have 
proposed a model in which cold relativistic pulsar wind (electron 
and/or proton loaded) with the bulk Lorentz factor T ~ 10" -r lO"* 
provides a source of X-ray emitting electrons. Such relativistic 
wind penetrating into the Be star disk can also be responsible for 
the formation of jet-like outflow from the system. 
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